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Abstract 

A fitting method to reconstruct the momentum and direction of charged particles 
in slightly inhomogeneous magnetic fields is presented in detail. For magnetic fields 
of the order of 1 T and inhomogeneity gradients as large as 1 T/m the typical 
momentum bias due to the proposed approximations is of the order of few MeV, 
to be compared with scattering components of the order of 20 MeV or even larger. 
This method is currently being employed in the reconstruction programs of the 
AMS experiment. 



1 Introduction 



The next generation of particle physics experiments will be characterized by 
an unprecedented accuracy in the determination of the positions and momenta 
of very energetic charged particles. The principle of momentum measurement 
is, in all cases, the linear relation between the track curvature and the inverse 
of the momentum in a plane perpendicular to the magnetic field direction. 
High energy physics experiments have successfully employed this principle for 
many decades, producing many relevant physics results and discoveries. 

Track fitting in inhomogeneous magnetic fields involves the propagation of 
track parameters between consecutive detector layers. Typical approaches [1,2], 
use numerical methods of high order in order to integrate the equations of mo- 
tion from layer to layer. This report describes a simple alternative algorithm, 
currently being employed in the AMS experiment [3]. All propagation opera- 
tions are expressed in terms of path integrals, which are approximated with 
enough accuracy at initialization. The fitting step is also reduced to a linear 
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problem. The simplicity of the algorithm allows fast refitting of tracks, even 
beyond the reconstruction phase, i.e. at the level of final physics analyses. 



The report is organized as follows. The approximations that are the basis of the 
method and few formulae to estimate its accuracy are presented in section 2. 
The track fitting logic and a typical implementation are described in section 3. 
Section 4 discusses the inclusion of multiple scattering effects. The report is 
summarized in section 5. 



2 Basis of the method 



The trajectory of a particle with charge q in a static magnetic field B is 
governed by the equation: 



f=?(*x3) (1) 

where p and v are the momentum and velocity of the particle at a given 
position x and time t on the trajectory. An immediate conclusion is that 
v = \\v\\ and p = \\p\\ are constants of motion. Locally, the trajectory is a 
helix, with p sin 9b = q\\B\\R, where R is the (signed) radius of curvature in a 
plane transverse to B and 9 b the angle between v and B. Equation 1 can be 
rewritten in a different way: 



dx v 
di v 

where dl = v dt is the differential length traversed by the particle. More visu- 
ally, u is a unitary vector tangent to the trajectory at the point x. Integration 
of this expression between two consecutive layers of a tracker detector, denoted 
by the subscripts and 1, leads to: 



m 1= m + ^ [ h dl'(ux B) (/') (3) 

P Jlo V ' 



The first approximation in our method consists in computing the previous 
equation as follows: 
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(4) 



that is, computing the integral along the straight line connecting xq and x±. 
The approximation is exact in two ) when the magnetic field is homo- 

geneous, and b) in the infinite momentum limit. For the homogeneous field 
case, B{1) = B , and Equation 3 becomes: 



Ui=U + - 

V 



( dl'uU')} x B = u + - (fi - x ) x B 
Jio J V 



(5) 



where the last equality is obtained by introducing the definition of u: dx = udl. 
The result is identical to the one obtained using Equation 4. In the limit of 



very high momentum, case b), u 
1/p, leading trivially again to Equation 4. 
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up to relative corrections of order 



In order to estimate the accuracy of the approximation in a general case, two 
additional expressions are necessary. First, the expression of the magnetic field 
as a series expansion around l\/2 = (/o + h)/^'- 



B(l)=B + B' (Z-Z 1/2 ) + ... 



(6) 



Second, the estimate A(-u) of the difference between the true u vector and 
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Introducing the previous expressions in Equation 3 and comparing with Equa- 
tion 4 one obtains the following correction at first order: 
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where R is the approximate radius of the trajectory from Iq to l\. The effect 
has to be compared with the corresponding term in Equation 4, of order 
|||.E>o||(/i — ^o)- The difference translates into a relative momentum shift of 
order: 



Ap 
P 



h - l AB 0.3 (h - Z )[m] AB[T] 



12 R 



x 



I fin 



12 p[GeV] 



(9) 



where AB is the typical variation of the magnetic field between l and l\. For 
instance, inhomogeneities in the magnetic field of the order of 1 T/m imply 
uncertainties of order Ap m ±1.6 MeV, independent of the absolute value of 
the rigidity. This has to be compared with the typical contributions from mul- 
tiple scattering in silicon detectors. For AMS-02 [3], optimized in this respect, 
the expected momentum resolution at the lowest momenta suggests multiple 
scattering effects of order Ap > 20 MeV [3], safely beyond the accuracy of the 
approximation. 

The second approximation concerns the extrapolation of the position vector 
onto the adjacent plane. Integrating Equation 2 twice we obtain: 



q rh rx , 
xi = ^o ± uo(h ~ k) + ~ / dx dy(uxB)(y) 

p Jlo JO ' 

= x + u (h - k) + - j dy(h ~ y) (u x B) (y) 

V Jin 



(10) 



The p — > oo limit reads: 



%i =x + Uo\\xi - foil + - / dy(l 

P Jlo 
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The previous expression, which is linear in x , u and q/p, does not coincide 
in general with the exact solution for the homogeneous magnetic field case. 
Nevertheless, we will prove that it is precise enough for most cases of interest. 
Let us define a convenient orthonormal reference system by the unitary vectors 
ua, ub and uq'- 
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U A = 



u B ; 



u c = 



(f i - f ) - [«b(^i - ^o)] Mb 



(fi - f )|| sm6>B 



B 
\B\\ 

(xi - x ) x m b 
|(fi - f )|| sin0 B 



(12) 
(13) 
(14) 



Note that u B is the unitary vector in the direction of the magnetic field, and 
9b the angle between the vectors (x\ — xq) and B. In terms of these vectors, 
the trajectory in a homogeneous field corresponds to: 



u(l) = \a(1)ua + \b{1)u b + Ac(/)«c; 
sin 6 b 



Aa(0 = sin#B cos 
As(/) = cos 6b 
^c(0 — sin^ B sin 



sin 6 b 
R 
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(15) 
(16) 
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with 6b and l re f constants. Let us also write (l± — 1 ) as an expansion in powers 
of l/R: 



n j \ I,- - ii . Il^i -x \\ 3 sm 2 d B , 
{h - l ) = \\ Xl -x \\ + + • • • 



(20) 



Equation 10 for the homogeneous case can be then rewritten as: 
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Since the linear term in q/p gives a contribution of order 



(h-k,) 2 

2R ' 



(21) 



a naive 



calculation would suggest a relative shift in the momentum of order: 



Ap h-l Q 0.3 (Zi-Zo)[m] 5[T] 
p ~ 12 i? 12p[GeV] 



(22) 
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For a benchmark separation of 20 cm and a magnetic field of 1 T, we obtain 
a maximum possible shift of Ap pa ±5 MeV. In practice, the effect is even 
smaller, since the missing correction affects coordinates in directions less sen- 
sitive to bending (u and ua). For most experiments, the fitting procedure 
is based on the minimization of a function in which position measurements 
in bending and non-bending directions are almost decoupled. In this config- 
uration, the correction above will act in quadrature, effectively leading to a 
momentum shift of order: 




0.3 (/!-/o)[m] B[T] 
12 p[GeV] 



(23) 



The quoted shift is negligible, even for large magnetic fields, like those of 
LHC and future linear collider detectors. In fact, this conclusion is somehow 
equivalent to the one reached in Reference [4] in the context of homogeneous 
magnetic fields. There, only measurements along the direction of the impact 
parameter with respect to the track at each point (i.e. the sensitive "bending" 
direction) are considered. This assumption leads naturally to a linear problem 
in terms of the curvature parameter [4]. 

In summary, the following approximations are considered to be accurate enough 
for most practical cases: 



uipau + -\\xi — x \\ [ da-^- — x B(x + a[x 1 — x ]) (24) 
p Jo \\Xi — x || 

XiPax + Uo\\xi — x \\ 

+ -\\x*i — x \\ 2 / da (1 — a) — ^ — x B(x + a[x\ — x }) (25) 

V JO \\Xi — Xn 



3 Track fitting 



For simplicity, it is assumed that all tracker sensitive layers are parallel to the 
z direction and that uncorrelated position measurements are performed along 
the x and y directions. The extension to more elaborated geometrical config- 
urations is straightforward, since only simple rotations of the predictions are 
involved. An obvious example is that of a detector with a radial configuration. 
To deal with it, it is enough to substitute one of the y 2 terms in the expressions 
presented later by a sum of residues along the azimuthal direction. 

We consider a scenario in which z coordinates are known with infinite pre- 
cision, so they can be fixed to their nominal values. The inclusion of an ad- 
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ditional z-term is, nevertheless, a trivial extension to the proposed scheme. 
Multiple scattering effects will be discussed in the next section. 

We need to determine the position of the track at the first plane, x = 
(x ,y ,z ), the tangent vector at the first plane, u = (u 0x , u 0y , u 0z ), and the 
inverse of the rigidity, q/p. From Equation 25 we obtain, on the second plane: 



Q ~* 2 

x i S3 f o + uoho + -Pmho 
p 

where the following definitions have been introduced: 



(26) 
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]) 



(27) 
(28) 



The integrals Pjj-i and the lengths are stored in an initialization phase. 
They are determined from the measured positions and the magnetic field val- 
ues on the line segment defined by Xj-i and Xj. 

At the third layer the extrapolation is given by: 



Q — * 
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(29) 



where the following path integral definition has been introduced (according to 
Equation 24): 
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x j x j-i 
\ x .i — x j-i 



x B(xj-i + a[xj — Xj-i]) 



(30) 



For many cases of interest (AMS-02 for instance), a Simpson method with just 
a few points is enough to calculate Pjj-i and jjj-i with sufficient accuracy. 
In general, the extrapolation to layer % can be written as: 




7 



where the lengths hk must be interpreted as follows: 



hk y ' h, 



m=k+l 



(32) 



The five parameters (x , y , u 0x , u 0y , q/p) are finally obtained by minimization 
of the following chi-square: 



2 _ ( X i,meas ~ Xjjxp, y , U 0x , U 0y , q/p))' 

^ ~ 2-^i „2 



i=0 °x 



_l_ (y i < meas ~ Vifaoi UOi u 0x, Up y , q/p)) 



i=0 S 



where 

%i,meas and Di^meas 

are the measured positions on layer i of the tracker, 
and a x and a y are the tracker position resolutions in the sensitive directions. 
The minimization leads to a linear equation, which can be easily solved via 
matrix inversion. 

Using an even more simplified notation, the x 2 can be written in a more 
convenient form: 



JV-l 



,- — \^i,meas X/fc=l ^ikVk) ^l—y \Vi,meas ^2k=l^*ikPk) 

X 2 =T,~ ~2 -+T,- ~2 (34) 

i=0 U x i=0 u y 

where p^j = 1,5 defines the vector of parameters to be determined, p = 
(x ,yo,Uo x ,u y,q/p). The components of the matrices Fij and are: 



Fn = l 
F i2 = 
Fi3 = ho 
F iA = 



F i5 = (Pk,k-Jl,k-i + lk,k~i l k,k-ihk) 



(35) 



k=i 



and: 
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{Pk,k-Jk,k-1 + lk,k~l l k,k-lkk) 



(36) 



k=i 



with the upper indices x and y denoting the x and y components of the vector 
integrals p k ,k-i and ~f k ,k-i- 



4 Multiple scattering treatment 

It will be assumed that the amount of traversed material is reasonably small 
and that the momentum range of interest is such that energy losses can be 
safely neglected. In these conditions, multiple scattering between layers j — 1 
and j is taken into account by estimating the additional uncertainty induced 
on the director vector Uj. This uncertainty depends on: a) the amount of 
traversed material in radiation lengths, b) the particle momentum and c) its 
velocity. 

From layer j — 1 to layer j a particle is traversing the amount of material Xjj-i, 
measured in radiation lengths. The rms angular deviations in the xz and yz 
projections are equal. Denoting them by by Ajj-i, they can be approximately 
parametrized [5,6] as follows: 



where (3 is the velocity of the particle (in c units) and p its momentum ex- 
pressed in GeV. Note also that the Xjj-i thicknesses hide a dependence on the 
director vectors Uj. The previous expression, accurate at the few percent level 
in the range 0.003 < -Xjj-i < 0.01 [6], does not admit a Gaussian treatment, 
in the sense that the expected additive property as a function of the amount 
of material is not satisfied: A 2 (X + Y) ^ A 2 (X) + A 2 (F). 

It is convenient to work with Gaussian uncertainties in order to keep a x 2 
minimization scheme. In the case of a very small amount of traversed material 
a possible approach is to assume the previous formula to be correct for the total 
amount of traversed material and then distribute the remaining deviations in a 
linear way at any intermediate plane, i.e. such that the rms deviation is always 
proportional to \fX. If the total amount of material is X tot , the suggestion 




(37) 



implies: 
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0.0136 [1 + 0.038 \n(X tot )} q 



p[GeV] 



(38) 



The previous estimate is usually consistent with the quoted accuracy of Equa- 
tion [6]. For the AMS-02 silicon tracker, it overestimates the rms deviations 
at the intermediate planes by at most 4%. 

From the fitting point of view, multiple scattering just modifies the directions 
at the different layers as follows: 



u^u^NO MS) + J2?k,k-i 
k=i 



(39) 



where Uj(NO MS) denotes the calculation in the absence of multiple scatter- 
ing and tjj-i is a deviation that follows a Gaussian of mean zero and width 
~ (Ajj-i, Ajj-i, 0). The j dependence enters through the amount of accumu- 
lated radiation lengths between the exit of layer j — 1 and the exit of 
layer j. At the level of position measurements the modified trajectories read: 



Xj=Xj(NO MS) + Cfc,fc-1 Yl l m+l,m = ?k,k-lljk 
k=l \m=k ) fc=l 



(40) 



There are two possible options to include these additional sources of un- 
certainty in the \ 2 ■ The first one is to fit all these additional parameters 
(2(number of planes — 2)) with additional Gaussian constraints according to 
the expected widths. We will employ a second option, keeping the same num- 
ber of fitted parameters, but building new covariance matrices according to 
the Gaussian uncertainties A fc)fc _ 1 . In the absence of multiple scattering, the 
covariance matrices for the x and y projections , V£ and are given by: 



V, 



al 
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• • err. 



(41) 



In the presence of multiple scattering, we need to take into account all fully 
correlated sources via terms, leading to the matrices: 
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min(i,j) — 1 

m=l 
min(i,j)—l 



m=l 



Finally, the x 2 reads: 



N-l I 5 \ 

X ^ ' I %i,meas ^ ' FikVk J 
i,j'=0 V fc=l / 

+ Z! (vi,meas ~ Z G ifcPfe) (vj 
i,j=0 \ k=l / \ 



or, in matrix form: 



x 



j,meas 



^ ' FjmPm 
m=l 

5 

^ ] GjmPr, 



m=l 



(42) 



(43) 



X 2 = (Xr. 



Fpf V~ l (x. 



Fp 



+ {Vmeas ~ Gp) W~ (y meas ~ Gp ) 



(44) 



Formally, its minimization with respect to p leads to the solution: 



p 



F T V^ 1 F + G T W^ 1 G 



F^^x^as + G'W-'y,, 



(45) 



Even if the x 2 to be minimized seems formally linear in the parameters p, 
multiple scattering introduces a dependence on q/p and uq via the covariance 
matrices V and W . A convenient way to solve the problem is to minimize 
the x 2 following an iterative procedure. In a first step, the x 2 is minimized 
using the diagonal covariance matrices and Wf,. The minimization is then 
iterated several times, using the Vij and Wij matrices determined from the 
parameters of the previous step. The iterative procedure is rapidly convergent. 
It may be stopped either after a couple of iterations or when some convergence 
criteria are reached. If computing time is not an issue, a convenient choice is 
to stop when the difference in rigidity between two consecutive steps is smaller 
than the accuracy of the method (a few MeV). 



5 Summary 



We have presented a simple algorithm for track fitting of high energy particles 
traversing slightly inhomogeneous magnetic fields. The method is based on the 
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prior calculation of a few path integrals which depend just on the measured 
positions and a few values of the magnetic field. The minimization of a x 2 , 
which presents a linear dependence on the track parameters, leads to a simple 
solution of the problem. Multiple scattering is considered in a straightforward 
and user-controlled way. This is particularly important when potential detec- 
tor resolution problems have to be disentangled from trivial material budget 
effects. Compared to other methods, a few simple formulae (9, 22 and 23) allow 
for a fast estimate of the expected momentum uncertainties. These formulae 
use as inputs the average value of the magnetic field, the typical size of the 
field inhomogeneities and the distance between measuring layers. The uncer- 
tainties translate into a shift of the measured momentum which, for most cases 
of interest in present and future high energy experiments, are of the order of a 
few MeV, well below the uncertainties due to multiple scattering. The method 
discussed here is being employed in the context of the AMS experiment [3]. 
Thanks to its intrinsic simplicity, it is also being used for fast and reliable 
track fitting at the latest steps of data reconstruction and analysis. 
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